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Abstract 

We report the measurements of nuclear modification factors of the Z bosons, isolated 
photons and charged particles in ^Jsm^ = 2.76 TeV PbPb collisions with the CMS de- 
tector. The nuclear modification factors are constructed by dividing the PbPb pj spec- 
tra, normalized to the number of binary collisions, by the pp references. No modifi- 
cations are observed in isolated photon and Z boson production with respected to the 
pp references while large suppression is observed in the charged particles. 
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1 Introduction 

The hot and dense matter produced in heavy-ion collisions, often referred to as the quark- 
gluon plasma (QGP), can be studied in various ways. The first way is to measure the yields of 
colorless particles in heavy-ion (AA) collisions that are unmodified in order to understand the 
initial state of the collisions. The second approach is to compare measurements made in AA 
collisions to those in proton-proton (pp) and proton- (or deuteron-) nucleus collisions. At the 
Relativistic Heavy Ion Collider (RHIC) direct photons play the reference role [Ij. At the Large 
Hadron Collider (LHC) energies, not only photons but also Z bosons, decaying to two leptons, 
become available El |3l. Since, once produced, photons and leptons from Z decays traverse 
the produced hot and dense medium without interacting strongly, they provide a direct test of 
perturbative QCD (pQCD) and the nuclear parton densities |4|. High transverse energy (Ej) 
prompt photons in nucleus-nucleus collisions are produced directly from the hard scattering 
of two partons. In order to suppress the much larger background coming from the electromag- 
netic decays of neutral mesons (mostly tz^, ?/ — > 7 7) produced in the fragmentation of the hard 
scattered partons, isolation requirements on the reconstructed photon candidates are imposed. 
The application of isolation requirements suppresses also the fragmentation photon component 
and enhances significantly photons coming from the quark-gluon Compton process [51. 

The measurement of charged particle pj spectrum is motivated by lower energy results H- 
m from the Relativistic Heavy Ion Collider (RHIC), w^here high-px particle production was 
found to be strongly suppressed relative to expectations from an independent superposition 
of nucleon-nucleon collisions. This observation is typically expressed in terms of the nuclear 
modification factor: 

where N^^ and cr^^^ represent the yield in nucleus-nucleus collisions and the cross section in 
nucleon-nucleon collisions, respectively. The pseudorapidity is defined as ;; = — ln[tan(0/2)], 
with 9 being the polar angle of the charged particle with respect to the counterclockwise beam 
direction. The nuclear overlap function (Taa ) is the ratio of the number of binary nucleon- 
nucleon collisions (NcoU ) calculated from a Glauber model of the nuclear collision geome- 
try [10] and the inelastic nucleon-nucleon cross section a^ = (64 ± 5) mb at i/s = 2.76 
TeV [|llil . In the absence of nuclear effects on the PbPb px spectrum, the factor Raa is unity 
by construction. In particular, together with recent studies of jet quenching and fragmentation 
properties HH [I3l, the modification of the px spectrum compared to nucleon-nucleon colli- 
sions at the same energy can shed light on the detailed mechanism by which hard partons lose 
energy traversing the medium [,14] . 

This report presents the measurements of nuclear modification factors of the isolated photons, 
Z bosons and charged particles as a function of event centrality for PbPb collisions collected 
by the Compact Muon Solenoid (CMS) experiment at a center-of-mass energy of 2.76 TeV per 
nucleon pair. 

2 Experimental methods 

Final-state particles in the selected collision events are reconstructed in the CMS detector. A 
detailed description of CMS can be found elsew^here |1T6| . A minimum-bias (MinBias) event 
sample is collected using coincidences between trigger signals from the +z and — z sides of 
either the Beam Scintillator Counters (BSC) or the Hadronic Forward Calorimeter (HF, covering 
2.9 < \^\ < 5.2) . The minimum bias trigger and event selection efficiency is 97 ± 3%. I15| . 



2 4 Results 

The events are also selected by the two-level trigger of CMS. The di-muon trigger requires 
two muon candidates in the muon detectors at the first hardware level (LI). At the software- 
based higher-level (HLT), two reconstructed tracks in the muon detectors are required, each 
with a Pj of at least 3 GeV/c. The photon trigger requires a LI electromagnetic cluster with 
Ej > 5 GeV and a HLT photon with Ej > 15 GeV. In order to extend the statistical reach of the 
charged particle pj spectra, data recorded by single-jet triggers with uncorrected transverse 
energy thresholds of Ej = 35 GeV and 50 GeV are included in the analysis. The jet-energy 
thresholds in the trigger are applied after subtracting the underlying event energy but without 
correcting for calorimeter response. In order to select a pure sample of inelastic hadronic col- 
lisions for analysis, the contamination from ultra-peripheral collisions and non-collision beam 
background are removed following the prescription of Ref. 1,15 J . 

3 Centrality determination 

For analysis of Fx)Pb events it is important to determine the overlap or impact parameter of 
the two colliding nuclei, usually called "centrality". Centrality is determined with the MinBias 
sample using the total sum of energy signals from both HF, where the events that deposit the 
top 10% most energy of the total interaction cross section are called the 0-10% most central 
(small impact paramter). Using Glauber model simulations, the intervals can be correlated 
with more meaningful physics quantities. Details of the centrality determination and the model 
calculations are described in Ref. [|15j . The two most commonly used centrality parameters 
are Npart/ the total number of nucleons in the two Pb nuclei (with mass number 208) which 
experienced at least one collision, and NcoU/ the total number of nucleon-nucleon collisions. 

4 Results 

4.1 Z boson 



The yield of Z — > }i^ji^ decays per MB event is defined as dN/(ii/(|i/| < 2.0) = Nz/{oceNMB^y)> 
where Nz is the number of dimuons counted in the mass window of 60-120 GeV/c^, Nmb = 
55 X 10^ is the number of corresponding MB events, corrected for trigger efficiency, a. and e are 
the acceptance and overall efficiency, and Ay = 4.0 is the rapidity bin width. The full circles in 
Fig. IT] (a) show the centrality dependence of the Z yield divided by Tab, ^vhile the open square 
is for MB events. The variable used on the abscissa is the average number of participating nu- 
cleons Npart corresponding to the selected centrality intervals, computed in the same Glauber 
model. No centrality dependence of the binary-scaled Z yields is observed in data [11 7) . A 
similar result was recently published by the ATLAS collaboration [18]. 

The normalized yields {dN/dy)/TAB are compared to various calculations [IT9l - l24l . Only a 
marginal centrality dependence is predicted: the inhomogeneous (i.e. depending on the radial 
position in nuclei) shadowing is predicted to have negligible impact II25II and the energy-loss 
prediction drops by 3% from peripheral to central collisions [22J. Figures [l] (b) and (c) show 
the differential yields, dN /dy and d^N /dydpj, as a function of the Z boson y and pj. The 
differential yields are compared to the same theoretical calculations as used for the central- 
ity distribution (when available) multiplied by the minimum bias Tab value. In all bins, no 
significant deviations from binary-collision scaling are observed. 

Nuclear modification factors, Raa = dN / [Tab x dcTpp), are computed from the AA measured 
yields dN, the nuclear overlap function Tab, and the pp— > Z cross sections dcTpp given by the 
POWHEG calculation (solid lines on Fig. [ll e.g. ddpp/dy = 59.6 pb in \y\ < 2.0). It is found that 
the Raa = 1-20 ± 0.29(stat.) ± 0.16(systO in 0-10% central PbPb collisions which is consistent 



4.2 Isolated photons 
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Figure 1: The yields of Z — > /^/y per event: a) dN/dy divided by the expected nuclear overlap 
function Tab and as a function of event centrality parameterized as the number of participating 
nucleons Npart/ b) dN/dy versus the Z boson y, c) d^N/dydpj versus the Z boson pj. Data 
points are located horizontally at average values measured within a given bin. Vertical lines 
(bands) correspond to statistical (systematic) uncertainties. Theoretical predictions are com- 
puted within the same bins as the data and are described in the text. 

with unity. 

4.2 Isolated photons 

The corrected photon yield with \f]\ < 1.44per MinBias event is defined as dN/dE^ = N'^/(e x 
Nmb X AEj), where Nub = 5.19 x 10'' is the number of MinBias events corresponding to the 
number of events recorded in the photon-triggered sample and e is the efficiency of the photon 
identification. The systematic uncertainty of the photon yield dN /dEj is dominated by the un- 
certainty of the background estimation which is in the level of 11-31% and the total systematic 
uncertainties are 21-37%. 

The nuclear modification factor Raa = dN'^ / {T^a ^ (^pp)i is computed from the PbPb mea- 
sured yield dW^, the inclusive isolated photon cross-section a^p given by the JETPHOX calcu- 
lation Il26ll . The uncertainty of Taa is included in the Raa systematic uncertainty. Figure |2] (left) 
shows the Raa as a function of the photon Ej in the 0-10% central collisions. The results are 
found to be compatible with unity within the quoted uncertainties. Figure [2] (right) shows the 
JETPHOX predictions for the same ratio using different nPDFs. 

Figure Islshows the Raa as a function of Npart in order to investigate the centrality dependence 
of the nuclear modification effect. No significant centrality dependence is observed in data. 

4.3 Charged particles 

The nuclear modification factor Raa is constructed according to Eq. ^ by dividing the PbPb re- 
sults by the scaled pp reference obtained from an interpolation [27J. Raa is presented as a func- 
tion of transverse momentum in Fig.|4]for each of the six centrality bins. In the most-peripheral 
events (70-90%), a moderate suppression of about two (Raa = 0.5) is observed at low px with 
Raa rising gently w^ith increasing transverse momentum. The suppression is increasingly pro- 
nounced in more-central collisions, as expected from the longer average path lengths traversed 
by hard-scattered partons as they lose energy via jet quenching. Raa reaches a minimum value 
of 0.13 around 6-7 GeV/c in the 0-5% bin. At higher pj, the value of Raa rises and levels off 
above 40 GeV/c at a value of approximately 0.5. A rising Raa niay simply reflect the flattening 
of the unquenched nucleon-nucleon spectrum at high pj, although the magnitude of the rise is 
not strongly constrained by theory (see Fig.lsll. 
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Figure 2: (Left Panel) Nuclear modification factor Raa as a function of the photon £7- 
measured in 0-10% central PbPb collisions over the T^A-scaled pp JETPHOX prediction at 
a/snn = 2.76 TeV. The vertical lines are the statistical uncertainty. The systematic uncertainties 
of the measurements, including the uncertainty of T^A/ are shown as yellow filled areas. The 
scale uncertainty of the pp calculation is shown as red lines. (Right Panel) The ratio of PbPb 
and pp JETPHOX predictions as described in the text is shown as the blue line. The uncertainty 
in this ratio due to the EPS09 PDF is shown as the blue dashed lines. Comparison between the 
JETPHOX predictions with alternative nPDFs is also shown. 



The evolution of the nuclear modification factor w^ith center-of-mass energy from SPS 
to RHIC ||30ll3ll to the LHC 133 is presented in Fig. |5) Note that results are shown for both 
charged hadrons and neutral pions, the latter being somewhat more suppressed at low pj (e.g., 
parton recombination can produce an excess of protons [33 1). At RHIC, the values of Raa from 
the two different measurements are seen to converge above roughly 8 GeV/c 11301 |3T'|. At low 
Pt/ the spectrum of charged particles is significantly more suppressed at the LHC than at RHIC, 
although the minimum values and the pj at which they occur are quite similar. It remains to 
be seen w^hether the similarity between the tz^ suppression observed by PHENIX at ^/s^ = 
200GeV [30] and the charged particle suppression presented in this paper at ^/s^ = 2.76TeV 
is more than a coincidence. 

The CMS measurement of Raa presented in this paper in the 0-5% centrality interval can be 
compared to the published ALICE [32j result over the pj range measured by ALICE. Note 
that the CMS pp reference [27J, derived from an interpolation that includes CDF and CMS 
measurements from ^/s = 0.63-7TeV, is roughly 5-15% higher than the ALICE pp reference 
quoted in their paper II32I . The tw^o results are in agreement within their respective statistical 
and systematic uncertainties. 

The high-px measurement of Raa from this analysis, up to pj = 100 GeV/c, is compared to a 
number of model predictions, both for the LHC design energy of a/s^ = 5.5TeV (PQM ["341 
and GLV [|35l[36l) and for the actual 2010 collision energy of ,/s;^ = 2.76TeV (ASW, YaJEM, 
and an elastic scattering energy-loss model with parameterized escape probability Il37li ). While 
most models predict the generally rising behavior that is observed in the data at high pj, the 
magnitude of the predicted slope varies greatly depending on the details of the jet quenching 
implementation. The new CMS measurement should help constrain the quenching parameters 
used in these models and further the understanding of the energy-loss mechanism. 
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Figure 3: The measured nuclear modification factor Raa as a function of Npart for the five 
different photon transverse energy intervals. The vertical lines are the statistical uncertainty. 
The systematic uncertainties of the measurements, including the uncertainty of T^A/ are shown 
as yellow filled areas. The scale uncertainty of the pp calculation is shown as red lines. 

5 Summary 



In conclusion, the Z boson and isolated photon yields in PbPb collisions at ^Snn = 2.76 TeV 
have been measured as a function of centrality Within uncertainties, no modification is ob- 
served with respect to theoretical next-to-leading order pQCD proton-proton cross sections 
scaled by the number of elementary nucleon-nucleon collisions. This measurement confirms 
the validity of the Glauber scaling for perturbative cross sections in nucleus-nucleus collisions 
at the LHC. 

A dramatic suppression of the charged particle spectrum has been observed for the most- 
central PbPb events, as was previously seen from the lower collision energies at RHIC. The 
result is consistent with the values of Raa measured by the ALICE experiment within experi- 
mental uncertainties. The rise in Raa from pj = 6-40GeV/c is reproduced in the predictions of 
a number of models, although the magnitude of the rise varies substantially between different 
predictions. Together with measurements of inclusive jet spectra, fragmentation functions, and 
energy balance, this measurement of the nuclear modification factor as a function of pj and 
collision centrality will help elucidate the mechanism of jet quenching and the properties of 
the medium produced in heavy ion collisions. 
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px-dependent systematic uncertainties on the Raa measurements. An additional systematic 
uncertainty from the normalization of Taa , common to all points, is shown as the shaded band 
around unity in each plot. 
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Figure 5: Measurements of the nuclear modification factor Raa in central heavy ion collisions 
at three different center-of-mass energies as a function of px for neutral pions and charged 
hadrons I.28.-31.I , compared to several theoretical predictions Il34l - l37il . The error bars around 
the points are the statistical uncertainties, and the yellow boxes around the CMS points are the 
systematic uncertainties. The bands for several of the theoretical calculations represent their 
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